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Abstract — We consider a generalization of Riemannian geometry that naturally 
arises in the framework of control theory. Let X and Y be two smooth vector fields 
on a two-dimensional manifold M. If X and Y are everywhere linearly independent, 
then they define a classical Riemannian metric on M (the metric for which they are 
orthonormal) and they give to M the structure of metric space. If X and Y become 
linearly dependent somewhere on M , then the corresponding Riemannian metric has 
singularities, but under generic conditions the metric structure is still well defined. 
Metric structures that can be defined locally in this way are called almost- Riemannian 
structures. They are special cases of rank-varying sub- Riemannian structures, which 
are naturally defined in terms of submodules of the space of smooth vector fields on 
M. Almost-Riemannian structures show interesting phenomena, in particular for what 
concerns the relation between curvature, presence of conjugate points, and topology of 
the manifold. The main result of the paper is a generalization to almost- Riemannian 
structures of the Gauss-Bonnet formula. 
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1 Introduction 



Let M be a two-dimensional smooth manifold and consider a pair of smooth vector fields X and Y on M. If 
the pair X, Y is Lie bracket generating, i.e., if span{X (q) , Y(q), [X, Y](q), [X, [X, Y]](q), . . .} is full-dimensional 
at every q G M , then the control system 

q = uX(q)+vY(q), U 2 + V 2 < 1 , q G M , (1) 

is completely controllable and the minimum-time function defines a continuous distance d on M. When X and 
y are everywhere linear independent (the only possibility for this to happen is that M is parallelizable), such 
distance is Riemannian and it corresponds to the metric for which (X, Y) is an orthonormal moving frame. 
The idea is to study the geometry obtained starting from a pair of vector fields which may become collincar. 
Under generic hypotheses, the set Z (called singular locus) of points of M at which X and Y are parallel is a 
one-dimensional embedded submanifold of M (possibly disconnected). 

Metric structures that can be defined locally by a pair of vector fields (X, Y) through (1) are called almost- 
Ricmannian structures. 

An almost-Riemannian structure can equivalently be seen as a locally finitely generated Lie bracket gen- 
erating C°°(M)-submodule A of Vec(M), the space of smooth vector fields on M, endowed with a bilinear, 
symmetric map G : A x A — > C°°(M) which is positive definite (in a suitable sense). A pair of vector fields X, 
Y in A is said to be orthonormal on some open set f2 if G(X, Y)(q) = and G(X, X)(q) = G(Y, Y)(q) = 1 for 
every ggft, 

An almost-Riemannian structure is said to be orientable if there exists a volume form, i.e., a bilinear, skew- 
symmetric, non-degenerate form ui : A x A — > C°°(M). In this case it is possible to normalize ui in such that 
\oj(X, Y)\ = 1 on fl for every open subset f2 of M and every local orthonormal frame (X, Y) on fl. 

It is interesting to notice that it is possible to define non-orientable almost-Riemannian structures on 
orientable manifolds and orientable almost-Riemannian structures on non-orientable manifolds. 

We say that an almost-Riemannian structure is trivializable if A is globally generated by a pair of vector 
fields defined on M. Trivializable almost-Riemannian structures are always orientable. 

The singular locus Z can be defined on M as the set where the linear subspace A(q) = {V(q) | V <G A} of 
T q M is not full-rank. An almost-Riemannian structure is Riemannian if and only if Z = 0. 

A famous example of genuine almost-Riemannian structure is provided by the Grushin plane, which is the 
almost-Riemannian structure on M 2 for which the vector fields X(x, y) = (1, 0) and Y(x, y) = (0, x) form a pair 
of orthonormal generators. (See Section 3.2 and [3, 4]. The model was originally introduced in the context of 
hypoclliptic operator theory [8, 9, 10, 11].) Notice that the singular locus is indeed nonempty, being equal to 
the y-axis. Another example of (trivializable) almost-Riemannian structure has appeared in problems of control 
of quantum mechanical systems (sec [5, 6]). In this case M = S 2 represents a suitable state space reduction of 
a three-level quantum system while the orthonormal generators X and Y are two infinitesimal rotations along 
two orthogonal axes, modeling the action on the system of two lasers in the rotating wave approximation (see 
Section 5.5). 

Such examples, and the naturalness of the construction leading to the definition of almost-Riemannian 
structure, motivate the study of general properties of such geometry, which exhibits many interesting features. 
One can check, for instance, that even in the case where the Gaussian curvature is everywhere negative (where 
it is defined, i.e., on M \ Z) geodesies may have conjugate points. For this reason it seems interesting to 
analyze the relations between the curvature, the presence of conjugate points, and the topology of the manifold 
(see also [1]). After providing a characterization of generic almost-Riemannian structures by means of local 
normal forms, in this paper we start this program by proving a generalization of the Gauss-Bonnet formula. 
Let M be compact and oriented, and endow it with an orientable almost-Riemannian structure. Denote by 
K ; M \ Z — ► ffi the Gaussian curvature. The first difficulty in order to extend the Gauss-Bonnet formula is 
to give a meaning to J M K dA, the integral of K on M with respect to the Riemannian density dA induced 
by the Riemannian metric on M \ Z. In the two examples cited above one can check, for instance, that, as q 
approaches Z, dA diverges, while K(q), which is everywhere negative, tends to — oo. 

The idea is to replace K dA with a signed version of it. A natural choice is K dA s , where dA s is a volume 
form intrinsically associated with the almost-Riemannian structure on M \ Z. 
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Our goal is to prove the existence and to assign a value to the limit 

lim f K(q)dA s , (2) 

e\0 J{ q& M\d(q,Z)>E} 

where d(-, •) is the distance globally defined by the almost-Riemannian structure on M. 

The goal will be attained under the following additional assumption. Generically the singular locus Z is 
smooth and A(q) is one-dimensional at every point of Z. We say that q € Z is a tangency point if A(q) is 
tangent to Z. Under generic assumptions, Z contains only a discrete set of tangency points. The hypothesis 
under which the main results of the paper are obtained is that Z contains no such point. 

Define M + (respectively, M~) as the subset of M \ Z on which the orientation defined by dA s coincides 
with (respectively, is opposite to) that of M. If M has no tangency point, then the limit (2) turns out to exist 
and is equal to 2ir(x(M + ) — x(M~)), where \ denotes the Euler characteristic. 

When the almost-Riemannian structure is trivializable, we have that x{M + ) = x(M~) and thus the limit 
(2) is equal to zero. Once applied to the special subclass of Riemannian structures, such result simply states 
that the integral of the curvature of a parallelizable compact oriented surface (i.e., the torus) is equal to zero. 
In a sense, in the standard Riemannian construction the topology of the surface gives a constraint on the total 
curvature through the Gauss-Bonnet formula, whereas for an almost-Riemannian structure induced by a single 
pair of vector fields the total curvature is equal to zero and the topology of the manifold constrains the metric 
to be singular on a suitable set. 

It is interesting to notice that every oriented compact surface can be endowed with a trivializable almost- 
Riemannian structure satisfying the requirement that there are no tangency points. 

The paper is organized as follows. In Section 2, we introduce two equivalent definitions of rank-varying 
sub-Riemannian structure on a manifold of any dimension, first by using the language of moduli and then by 
identifying it with an atlas of orthonormal frames. Rank- varying sub-Riemannian structures have already been 
studied, from a different perspective, in [12, 15]. 

A notion of orientability for rank-varying sub-Riemannian structure is then introduced. Almost-Riemannian 
structures arc defined as rank-varying sub-Riemannian structures of maximal rank. 

Starting from Section 3 we focus on the case of almost-Riemannian structures on two-dimensional manifolds. 
Geodesies associated with such structures are characterized in Section 3.1 using the Pontryagin Maximum 
Principle. In Section 3.2 we study the Grushin plane, for which we compute the cut and the conjugate loci. 

In Section 4 we provide local normal forms for generic almost-Riemannian structures, which are used in 
Section 5 to prove a generalization of the Gauss-Bonnet formula to almost-Riemannian structures without 
tangency points. The formula is then specialized to the case of trivializable almost-Riemannian structures. In 
Section 5.5 we show that every compact orientable two dimensional manifold admits a trivializable almost- 
Riemannian structure with no tangency points. 

2 Rank-varying distributions and sub-Riemannian structures 

Let M be a n-dimensional smooth manifold. Recall that Vec(M), the set of smooth vector fields on M, is 
naturally endowed with the structure of C°°(M)-module. Given an open subset fl of M, a submodule A of 
Vec(M) is said to be generated on f2 by the vector fields {V\, . . . ,V m } if every V £ A can be written as 
V = a\V\ H h a m V m on fl where ai,...,a m belong to C°°(M), 

Definition 1 A (n, k) -rank-varying distribution is a pair (M, A) where M is a n- dimensional smooth manifold, 
A is a submodule of Vec(M), and k < n is such that for every q G M and every small enough neighborhood 
n q of q, the restriction to SI 9 of A is generated by k vector fields and cannot be generated by less than k vector 
fields. 

From now on the expression A(q) will denote the linear subspace {V(q) | V € A} C T q M. When the 
dimension of A(q) is independent of q, we recover the standard definition of distribution as a smooth field of 
linear subspaccs of T q M. Notice that A cannot be identified with the map q i— > A(q). Indeed, it can happen 
that two different moduli Ai and A2 are such that Ai(q) = A2(q) for every q € M. Take for instance M = ffi 
and Ai, A 2 generated, respectively, by Fi(x) — x, F 2 {x) — x 2 . 
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Denote by Lie(A) the smallest Lie subalgebra of Vcc(M) containing A and let Lie g (A) = {V(q) | V £ 
Lic(A)} for every q £ M. We say that (M, A) satisfies the Lie bracket generating condition if Lic g (A) = T q M 
for every q £ M. We also introduce the flag of a rank-varying distribution (M, A) as the sequence of submodules 
A = A C Ai C . . . C A m C • • • defined through the recursive formula 

A fc+ i = A fc + [A,Afc]. (3) 

As above, we let A m (q) = {V(q) \ V G A m }. 

In order to provide an example of rank- varying distribution, let us introduce the Grushin distribution. Take 
as M the plane M 2 and let A be generated by the vector fields Fi(x,y) ~ (1,0) and F2(x,y) = (0,x). Then 
(K 2 , A) is a (2, 2)-rank-varying distribution. Notice that A(g) is equal to K x {0} when q is on the y-axis and 
to R 2 elsewhere. The Grushin distribution is Lie bracket generating since ^(q) = R 2 for every q £M. 2 . 

Crucial in what follows is the notion of generic (n, fc)-rank-varying distribution. Denote by W the C 2 - 
Whitney topology defined on Vec(M) and by (Vec(M), W) k the product of k copies of Vec(M) endowed with 
the corresponding product topology. We recall that if M is compact (as it is the case in most of what follows), 
then W is the standard C 2 topology. 

Definition 2 A property (P) defined for (n, fc) -rank-varying distributions is said to be generic if there exists 
an open and dense subset O of (Vec(M),W) fc such that (P) holds for every (n, fc)- rank- varying distribution 
which is generated by elements of O. 

E.g., gencrically, a (n, fc)-rank-varying distribution is Lie bracket generating. 
2.1 Orientable rank- varying distributions 

Let (M, A) be a (n, fc)-rank-varying distribution. A fc-form on (M, A) is a multilinear skew-symmetric map 

u : A x ■ ■ ■ x A -> C°°{M). 

k times 

We say that a fc-form oj is a volume form if, for every q G M, there exist k vector fields Fi, . . . ,Fk G A such 
that a}{F 1 ,...,F k )(q) ^0. 

Definition 3 We say that a (n, k) -rank-varying distribution (M, A) is orientable if it admits a volume form, 
otherwise we say that (M, A) is non- orientable. 

Notice that a rank-varying distribution can be orientable even if M is a non-orientable manifold (see example 
below). However, the distribution (M, Vec(M)) is orientable if and only if M is. 

Remark 4 Thanks to its multilinearity, a volume form is completely characterized by its action on the genera- 
tors. Given a trivializable rank- varying distribution A and a global system of generators F\, . . . , Fk, the equality 
. . . , Fk) = 1 uniquely defines a volume form on A. Hence every trivializable rank- varying distribution is 
orientable. 

Remark 5 Let ubca fc-form on a (n, fc)-rank- varying distribution (M, A). Then u> acts as a tensor on the open 
subset of M made of points q such that the dimension of A(g) is equal to fc, i.e., for every V\, . . . , Vk G A the 
value of u>(Vi, . . . , Vk)(q) depends only on Vi(q), . . . , Vk(q). Indeed, let {Fi, . . . , Fk} be a local system of gener- 
ators of A on a neighborhood ft of q and take ay £ C°°(M), 1 < i, j < fc, such that V{ = ^2j=i a ij^i 011 ^- The 
multilinearity and skew-symmetricity of lu implie that w(Vi, . . . , Vk)(q) = det(aij(q))u)(Fi, . . . , Fk)(q). There- 
fore, ui(Vi, . . . , Vk){q) depends on V\, . . . , Vk only through the matrix {aij(q)), which is uniquely determined by 
V 1 (q),...,V k (q). 

Let us present some example of orientable and non-orientable rank- varying distributions. All these examples 
are (2, 2)-rank-varying distributions, since this is our main case of interest in the following. 
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The Grushin distribution. Let M = K 2 and recall that A is generated by the vector fields F\(x, y) = (1, 0) 
and F2(x,y) = (0,x). The distribution A is orientable. A volume form can be defined by 

u(Vi,V 2 )(x,y) = lim Ijdx A dy (Vi(x' , y), V 2 (x' , y)), 

x'—>x X 

for every pair of vector fields Vi, V2 belonging to A. Equivalcntly we could have defined uj on the generators 
Fi, F 2 by requiring that ui(Fi, F 2 )(x,y) = 1. 
Notice that u> is not a tensor on the y-axis. 

A non-orientable rank- varying distribution on the torus. Let M be the two-dimensional torus [— tt, tt] x 
[— tt, tt] with the standard identifications. Consider the open covering of M given by 

Q 1 = (-it/2, tt/2) x [-tt, tt], n 2 = ([-tt, -7r/4) U (tt/4, tt]) X [—TT, TT]. 

Let A be generated by the vector fields 

2^ = (1,0), Fl = (0,sinx), onO 1 
2? =(1,0), F 2 ^ (0,1), onJJ 2 . 

This rank-varying distribution is non-orientable. In fact a volume form ui should acts on the local generators 
as 

u(Fl,Fl)(q) = fi(q) for every q G Q 1 , u(F 2 , F 2 )(q) - f 2 (q) for every q G O 2 , 
where fi and f 2 are two never- vanishing smooth functions. On O 1 n O 2 we would have 
fi(x,y)=u(FlF})(x,y)=sm(x)u>(FlF*)(x,y)=sm(x)f 2 (x,y), 
contradicting the constant-sign assumption on fi, f 2 . As a consequence A is not trivializable. 

An orientable rank-varying distribution on the Klein bottle. Let M be the Klein bottle seen as 
the square [— tt,tt] x [— tt,tt] with the identifications (x, — tt) ~ (x,tt), (— tt, y) ~ (tt, — y). Consider the vector 
fields 

Fi(x,y) - (1,0), F 2 (x,y) = (0, sin(2z)), 
which are well defined on M. The distribution generated by Fi and F 2 is orientable since it is trivializable. 



2.2 Rank-varying sub-Riemannian structures 

In this section we see how to introduce a a smoothly- varying Riemannian structure on every subspace A (q). 

Definition 6 A (n, k) -rank-varying sub-Riemannian structure is a triple S = (M, A, G), where (M, A) is a Lie 
bracket generating (n, k) -rank-varying distribution and G : A x A — > C°°(M) is a symmetric, positive definite 
bilinear map, i.e., a map such that for every V, W G A and f G C°°(M) we have 

G(V,W) = G(W,V), 
G(fV, W) = G(V, fW) = fG(V, W), 
G(V, V)(q) > 0, for every q G M, 
G(V, V)(q) = implies that V(q) = 0. 

A (n,n) -rank-varying sub-Riemannian structure is called a n-dimensional almost- Riemannian structure. 
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Let iS = (M, A, G) be a (n, fc)-rank-varying sub-Ricmannian structure. Reasoning as in Remark 5, we get that 
G is a tensor at the points q where dim(A(g)) = k. Although this is not necessarily the case everywhere on 
M, we can define, for every q G M, a quadratic form G q on A(g) through 

G q (v, v) = m£{G(V, V)(q) \ V(q) = v, V £ A}. 

For every q G M, it is possible to find a neighborhood Q q of g and an orthonormal frame on fi g , i.e., a set 
of fc vector fields X\,..., G A such that G(Xi, Xj) = Sij on f2 9 . One easily proves that orthonormal frames 
are local generators in Q q . 

If S is orientable then a volume form uj can be chosen in such a way that |w(Ai, . . . , X/.)\ = 1 on every local 
orthonormal frame. 

Let S = (M, A, G) be a (n, fc)-rank-varying sub-Riemannian structure. A curve 7 : [0, T] — ► M is said to be 
admissible for S if it is Lipschitz continuous and jit) G A 7 ( t ) for almost every t G [0, T]. Given an admissible 
curve 7 : [0, T] — * M, the length of 7 is 

J( 7 ) = ^ ^/G 7(t) ( 7 (t),7W) 

The distance induced by S on M is the function 

d(qo,Qi) = inf{^( 7 ) | 7 (0) = q ,l(T) = q lt -y admissible}. (4) 

It is a standard fact that ^( 7 ) is invariant under reparameterization of the curve 7. Moreover, if an admissible 
curve 7 minimizes the so-called energy functional £'(7) = Jq G 1 ^ t j( A /(t), j(t)) dt with T fixed (and fixed initial 
and final point) then v = y/G 7 ( t - ) ( A f(t),j(t)) is constant and 7 is also a minimizer of l(-). On the other side a 
minimizer 7 of /(•) such that v is constant is a minimizer of E(-) with T = l(p/)/v. 

A geodesic for S is a curve 7 : [0,T] — ► M such that for every sufficiently small interval [ti,ta] C [0, T], 
7l[ti.t 2 ] i s a minimizer of E(-). A geodesic for which G 7 ( t )( 7 (t), 7 (t)) is (constantly) equal to one is said to be 
parameterized by arclength. 

The finiteness and the continuity of d(-, •) with respect to the topology of M are guaranteed by the Lie 
bracket generating assumption on the rank- varying sub-Ricmannian structure. The distance d(-, ■) gives to M 
the structure of metric space. The local existence of minimizing geodesies is a standard consequence of Filippov 
Theorem (see for instance [2]). When M is compact any two points of M are connected by a minimizing geodesic. 

A convenient way to deal with a rank-varying sub-Riemannian structure is to identify it with an atlas of 
local orthonormal frames. In the case of an orientable rank-varying sub-Ricmannian structure, one can impose 
that all orthonormal frames are coherently oriented. In this way one is led to the following equivalent definition. 

Definition 7 Let M be a n-dimensional smooth manifold, fix k £ N, and consider a family 

S = {(f^Af,...,At)W, 

where {il^j^g/ is an open covering of M and, for every /1 G /, {Aj 1 , . . . , X^} is a family of smooth vector fields 
defined on M , whose restriction to satisfies the Lie bracket generating condition. We assume moreover that 
for every (i£j and every open nonempty subset of f2 M , the submodule of Vec(f2) generated by X^, . . . ,X£ 
on fl cannot be generated by less than k vector fields. 

We say that S is a (n, k) -rank-varying sub-Riemannian structure if, for every \x,v G / and for every 
q G H Q" , there exists an orthogonal matrix R^' u (q) = (Rfj{q)) G 0(k) such that 

k 

x?(q) = Y,R?;?(q)Xj(v)' (5) 

We say that two rank-varying sub-Riemannian structures Si and S2 on M are equivalent if S\ U S2 is a rank- 
varying sub-Riemannian structure. Given an open subset £1 of M and a set of k vector fields (Xi, . . . ,X/.), we 
say that (f2, X\, . . . , Xk) is compatible with S if S U {(fi, X\, . . . , X^)} is equivalent to S. 
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If S is equivalent to a rank-varying sub-Riemannian structure of the form {{M, X\, . . . ,Xk)}, i.e., for which 
the cardinality of I is equal to one, we say that S is trivializable . 

If S admits an equivalent rank-varying sub-Riemannian structure such that each belongs to SO(k), 

we say that S is orientable. 

Notice that R fl ' u (q) is uniquely denned by equation (5) and, moreover, is smooth as a function of q. In the 
following, when dealing with an orientable rank-varying sub-Riemannian structure, we always assume that the 
atlas of local orthonormal frames is positive oriented, i.e., such that each R' 1 ' 1 ' belongs to SO(k). For such an 
atlas, a volume form ui can be chosen such that u>(X^, . . . , X£) = 1 on f^, for every fi E I. 

In terms of Definition 6, A is the module that is locally (in f2 M ) generated by X^, . . . ,X£. Moreover 

A(q) = sv<m{X?(q),...,X%(q)}, 

{k k ^ 

Y,<x 2 i \v = Y,a i X?(q)\, 
i=\ i=\ ) 

for every q G Of 1 and v G A(g). 

Definition 8 A property (P) defined for (n, k) -rank-varying sub-Riemannian structures on M is said to be 
generic if there exists an open and dense subset O of (Vec(M), W) k such that (P) holds for every (n, ^-rank- 
varying sub-Riemannian structure admitting an atlas of local orthonormal frames whose elements belong to 
O. 

Given a (n, fc)-rank-varying sub-Riemannian structure iS, the problem of finding a curve minimizing the 
energy between two fixed points qo, q\ G M is naturally formulated as the optimal control problem 

k 

q = Y, Ui X?(q), u t eR, iiel{q)={iiel\qeW}, (6) 

2 = 1 

7 1 k 

^2 u i (*) dt min > ?(°) = Qo, q(T) = q\. (7) 

i=l 

Here /U, U\, . . . , are seen as controls and T is fixed. It is a standard fact that this optimal control problem is 
equivalent to the minimum time problem with controls u\, . . . , satisfying u\ + • • • + u\ < 1. 

Notice that if the rank-varying sub-Riemannian structure is trivializable, then the role of \x is empty and 
(6), (7) can be rewritten as a classical sub-Riemannian control problem 

q = y2 u i x i(q) , u i e R ) / y2 dt ~* min > 9(0) = qo, q(T) = q\. 

<=i Jo i=i 



3 Two-dimensional almost-Riemannian structures 

Henceforth the paper is focused on the special case of (2, 2)-rank-varying sub-Riemannian structures, i.e., two- 
dimensional almost-Riemannian structures (2-ARSs for short). In this case a local orthonormal frame on fi M 
is a pair of Lie bracket generating vector fields (X^, Y^). 

Given a 2-ARS S, we call singular locus the set ZcMof points q at which the dimension of A(g) is equal 
to one. Denote by g the restriction of the quadratic form G on M \ Z. By construction g is a Riemannian 
metric satisfying 

g(X»(q),X»(q)) = 1, g(X»{q),Y»(q)) = 0, g(Y»(q),Y»(q)) = 1, 

for every /i in I and every q G f2 M \ Z. Denote moreover by dA the Riemannian density associated with 
(M \Z,g), which coincides with \dX^ A dY^\ on il M \ Z, for every /16/. 
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Finally, one can define on M\Z the Gaussian curvature K associated with g, which is easily expressed in 
each open set f2 M \ Z through the formula (see for instance [2], equation (24.6)) 

K = -(o^) 2 - (Z^ 1 ) 2 + A^ - Y^a^, 

where , /3 M : fi^ \ Z — * R are (uniquely) defined by 

[X> t ,Y tl ] = afX* + /3^, 

and X^pv (respectively, Y^a^) denotes the Lie derivative of (3^ with respect to A M (respectively, of a? with 
respect to Y^). 

3.1 Minimizers, cut and conjugate loci 

A natural tool to look for geodesies in almost-Ricmannian geometry is to apply the necessary condition for 
optimality given by the Pontryagin Maximum Principle (see [13]). As a result we obtain the following proposi- 
tion. In view of later applications in the paper, we consider as initial and final conditions not only points, but 
submanifolds as well. 

Proposition 9 Define on T*M the Hamiltonian 

H(X,q) = ±((X,X»(q)) 2 + (X,Y^q)) 2 ), q&W, X&T*M. 
(Notice that H is well defined on the whole T*M, thanks to (5).) Consider the minimization problem 

qeA(q), [ G q(t) (q(t),q(t)) dt ^ min, q(0) = M in , q(T) = M fin , (8) 
Jo 

where M- m and Mfi n are two submanifolds of M and the final time T > is fixed. Then every solution of 

(8) is the projection on M of a trajectory (\(t),q(t)) of the Hamiltonian system associated with H satisfying 
A(0) J_ T m M in , A(T) J_ T 9(T) M fin; and H(X(t),q(t)) ? 0. 

Remark 10 The simple form of the statement above follows from the absence of abnormal minimizers, which 
follows from the Lie bracket generating assumption. As a consequence a curve is a geodesic if and only if it the 
projection of a normal extremal. 

Remark 11 Notice that H is constant along any given solution of the Hamiltonian system. Moreover, H = 1/2 
if and only if g(.) is parameterized by arclength. 

Fix q e M. For every A € T*M satisfying 

ff(A,g) = l/2 (9) 

and every t > define E(X,t) as the projection on M of the solution, evaluated at time t, of the Hamiltonian 
system associated with H, with initial condition A(0) = A and q(Q) = q. Notice that if q ^ Z then condition 

(9) defines an ellipse in T*M ; otherwise it identifies a pair of parallel straight lines. 

Definition 12 The conjugate locus from q is the set of critical values of the map E(X,t). For every A such 
that (9) holds, let t(X) be the first positive time, if it exists, for which the map (A, t) i— ► E(X,t) is singular at 
(A,£(A)). The first conjugate locus from q is the set {E(X,t(X)) | i(A) exists}. The cut locus from q is the set 
of points reached optimally by more than one geodesic, i.e., the set 

{q G M | 3 A 1; A 2 , t such that q' = E(X 1 ,t) = E(X 2 ,t), Ai ^ A 2 , and E(Xi, •), E(X 2 , ■) are optimal in [0,t]}. 

Remark 13 It is a standard fact that for every A satisfying (9), the set T(A) = {i > | the map (A, t) i— ► 
E(X,t) is singular at (A,t)} is a discrete set (see for instance [2]). 
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Figure 1: Geodesies and minimum time front (for t = 1) for the Grushin metric 



3.2 An example: the Grushin almost-Riemannian structure 

Consider again the Grushin distribution X(x,y) = (1,0), Y(x,y) = (0,x) on the plane R 2 . If we consider X 
and Y as an orthonormal frame, we get an almost-Riemannian structure. 

As already remarked, the singular locus coincides with the y-axis. Therefore, every trajectory crossing the 
y-axis does it horizontally. The Riemannian metric g associated with the Grushin metric on ffi 2 \ {(x, y) £ M. 2 | 
1^0} explodes when one is approaching the y-axis, 

g = dx 2 + \dy 2 . 

x z 

Also the curvature and the Riemannian density explode while approaching the y-axis, 

2 1 

K = -, dA=—dxdy. 

x* \x\ 

According to Proposition 9, geodesies are the projection on the (x, y)-plane of the solutions of the Hamil- 
tonian system corresponding to 

H = ±(\l + \ 2 y x 2 ). 

Fixing the initial condition x(0) = 0, j/(0) = 0, the normalization H = 1/2 implies that A x (0) = ±1. Taking 
a = X y (0) G K, the geodesies starting from the origin are 




Due to the symmetries of the problem, one can easily check that the time at which a geodesic (x a (t),y a (t)) 
loses optimality is t = 7r/|o|, for a ^ 0, and that (x a (t), y a (t)) belongs to the y-axis. The geodesies corresponding 
to a = are optimal for every positive time. As a consequence the cut locus from the origin is the set 
{(0,a) |aeR\{0}}. 
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In Figure 1A geodesies for some values of a are portrayed, while Figure IB illustrates the set of points 
reached in time t = 1. Notice that this set has a non-smooth boundary. In contrast with what would happen in 
Riemannian geometry, this is the case for every positive time, as it happens in constant-rank sub-Ricmannian 
geometry. However, this is a consequence of the fact that the initial condition belongs to Z. 

One can check that, even if the curvature is always negative where it is defined, a geodesic (x a (t),y a (t)), 
a 7^ 0, has its first conjugate point at time r/|a|, where r ~ 4.49 is the first positive root of the equation 
tan(r) = r. As a consequence the first conjugate locus is the parabola 

,.2 



X 2 1 1 



y = 



2 \ cos r sin r r 

One could ask whether the presence of conjugate points is the consequence of the particular initial point on 
the set Z. In fact this is not the case. Consider as initial condition the point x(0) = —1, y(0) = 0. Define, for 
every a S [0, 1], 



x + (t, a) 
x~ {t, a) 
y + (t,a) 
y~{t,a) 



— (a cos(ai)) + \J\ — a 2 sin(at) 



(a cos(ai)) — Vl — a? sin(ai) 



a 

\2 



-4 a Vl - a 2 + 2 a t + 4 a \/l - a 2 cos(a t) - sin(2 at) + 2 a 2 sin(2ai) 

2 



4a Vl - a 2 + 2 at - 4a Vl - a 2 cos(at) - sin(2at) + 2 a 2 sin(2 at) 



4 a 2 

Then every geodesic from the point (—1,0) belongs to one of the four families 



Gl 

G2 
G3 
G4 



{x + (t,a),y + (t,a)) } 
(x~(t,a),y~(t,a)), 
(x+(t,a),-y+(t,a)), 
{x-{t,a),-y + (t,a)). 



The geodesies in Gl and G2 lie in the half plane {y > 0}, while those in G3 and G4 lie in {y < 0}. 

Let us describe the cut locus from (—1,0). Consider first a geodesies (x + (t, a), y + (t, a)) belonging to the 
family Gl. One can check that ir/a is the first positive time at which (x + (t, a), y + (t, a)) intersects another 
geodesic, namely (x~ (t, a) , y~ (t, a)) , which belongs to G2. The situation is similar for the families G3 and G4. 
As a consequence the cut locus from (—1,0) is the set {(l,a) | a £ [tv/2, cc) U (— oo, — 7r/2]}. 

As above, one can also check that every geodesic (except those corresponding to a = 0) has a conjugate 
time (see Figure 2). In particular it is a conjugate time for the geodesies corresponding to a = 1. Notice 
that conjugate points appear on geodesies which have already crossed Z. (Before crossing Z a geodesic is 
Riemannian and lies in a Riemannian space with negative Gaussian curvature.) 



4 Normal forms for generic 2-ARSs 

The following proposition is a standard corollary of the transversality theorem. It formulates generic properties 
of a 2-ARS in terms of the flag of the distribution A (see equation (3)). 

Proposition 14 Let M be a two-dimensional smooth manifold. Generically, a 2-ARS S = {(f^, AT M , Y^)}^ 
on M satisfies the following properties: (i) Z is an embedded one- dimensional smooth submanifold of M ; (ii) 
The points q € AI at which ^2(0) is one- dimensional are isolated; (iii) A^(q) = T q M for every q £ M . 

Remark 15 Notice that properties (i), (ii), and (iii) are actually generic for every (2, 2)-rank- varying distri- 
bution, since they do not involve the metric structure. 
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Cut Locus 




Figure 2: 



As a consequence of Proposition 14. one can classify the local normal forms of a generic 2-ARS. 

Theorem 16 Generically for a 2-ARS S , for every point q £ M there exist a neighborhood U of q and a pair 
of vector fields (X,Y) on M such that (U,X,Y) is compatible withS and, up to a smooth change of coordinates 
defined on U , q = (0, 0) and (X, Y) has one of the forms 

(Fl) X(x,y) = (1,0), Y(x,y) = (0,e^), 

(F2) X(x,y) = (1,0), Y(x,y) = (0,xe^), 

(F3) X(x,y) = (1,0), Y(x,y) = (0,(y-x 2 ^(x))e^), 

where 4> and ip are smooth real-valued functions such that <p(0, y) = and t/j(0) =/= 0. 
Before proving Theorem 16 let us show the following lemma. 

Lemma 17 Let S be a 2-ARS and let W be a smooth embedded one- dimensional submanifold of M . Assume 
that W is transversal to the distribution A, i.e., such that A(q) + T q W = T q M for every q £ W. Then, for 
every q £ W there exists an open neighborhood U of q such that for every e > the set 

{q' £U\d(q',W)=e), 

is a smooth embedded one- dimensional submanifold of U . Moreover, there exists a pair of vector fields (X, Y) 
such that (U,X,Y) is compatible with S and, up to a smooth change of coordinates defined on U , q — (0,0) 
and W,X, Y have the form 

W = {(0,h) | ftGR}, 
X(x,y) = (1,0), 
Y(x,y) = (0,f(x,y)) 

where f(x,y) is a smooth function defined on U. 
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Proof. Consider a smooth regular parametrization a i— > w(a) of W. Let a i— > \o(a) G T^,^M be a smooth 
map satisfying H(X (a),w(a)) = 1/2 and A (a) _L T^^W. 

Let E(t, a) be the solution at time t of the Hamiltonian system given by the Pontryagin Maximum Principle 
with initial condition (q(0), A(0)) = (w(a), Ao(oj)) (see Proposition 9). Fix q G W and define a by q = w(a). In 
order to prove that E(t, a) is a local diffeomorphism around the point (0, a), let us show that the two vectors 

Vi = —{0,a) and u 2 = -^(0, a) 

are not parallel. On one hand, since v\ is equal to ^j(a), then it spans T q W . On the other hand, being H 
quadratic in A, 

OH 

(A (a),t> 2 ) = (A (a),— (A (a),?)> =2H{X (a),q) = 1. 

Thus V2 does not belong to the orthogonal to Ao(a), that is, to T q W . 

Therefore for a small enough neighborhood U of q we have that the set {q' € U \ d(q' , W) = e} is given by 
the intersection of U with the images of E(e, •) and E(—e, •). This proves the first part of the statement. To 
prove the second part, let us take (t, a) as a system of coordinates on U and define the vector field X by 

dE{t,a) 
X(t,a) = — gj— . 

Notice that, by construction, for every q' G U the vector X(q') belongs to A(g') and G q >(X(q'), X(q')) = 1. In 
the coordinates (t, a) we have X = (1, 0). Let Y be a vector field on U such that (X, Y) is compatible with S. 
We are left to prove that the first component of Y is identically equal to zero. Indeed, were this not the case, 
the norm of X would not be equal to one. 



Proof of Theorem 16. Let us start from the case in which A(q) = T q M . Let W be any one-dimensional 
submanifold passing through q. Lemma 17 provides us with a possible choice of orthonormal frame X = (1, 0), 
Y = (0, f(x,y)) in a neighborhood U of q. Since, without loss of generality, X and Y are everywhere linearly 
independent in U, then f(x,y) ^ for every (x,y) G U. By applying a smooth coordinate transformation of 
the type x — > x, y — > v{y) we get the new expressions X = (1,0), Y = (0, v'(y)f(x, y)). A normal form of type 
(Fl) is obtained by choosing v in such a way that i/(y)/(0, y) = 1. 

Let now q G Z and assume that ^(q) = T q M . Assume, moreover, that the generic condition (i) holds 
true. One can easily check that A(g) is transversal to the submanifold Z at q. Hence we can apply Lemma 
17 with W = Z. As a result we obtain a possible choice of orthonormal frame X = (1, 0), Y = (0, f(x, y)) in 
a neighborhood U of q. Since X and Y are linearly dependent on Z, which is identified with the y-axis, then 
f(0,y) = 0. The condition A 2 (q) = T q M implies that, by taking U small enough, d x f(0,y) ^ 0. Hence / 
admits a representation of the type f(x,y) = xe^ x,y \ with <f> smooth. Again, a change of coordinates x — ► x, 
y — » v{y) can be used in order to ensure that </>(0, y) = 0. The normal form (F2) is obtained. 

Let now q G Z be such that A2(q) = A(g). Assume that the generic conditions (i), (ii), (iii) are fulfilled. 
Let W be any one-dimensional submanifold passing through q and being transversal to Z. Using Lemma 17 
we can chose X — (1,0), Y = (0,f(x,y)) with / satisfying, by assumption, 9^/(0,0) = 0, 9^/(0,0) ^ 0. 
Let us identify Z with the graph of a smooth function y = T(x). Then f(x,y) can be written in the form 
(y — T(x))e^ x ^ with </> smooth. As above, without loss of generality </>(0, y) = 0. The conditions on / at (0, 0), 
moreover, justify the representation T(x) = x 2 ip(x), with ip smooth. 

Remark 18 Because of Remark 15, for a generic distribution A, every 2-ARS having A as corresponding 
distribution can be locally represented by one of the normal forms (Fl), (F2), (F3). 

Definition 19 Let S be a 2-ARS and assume that the generic conditions (i), (ii), (iii) of Proposition H hold 
true. A point q G M is said to be an ordinary point if A(q) = T q AI, hence, if S is locally described by (Fl). We 
call q a Grushin point if A(q) is one- dimensional and A.2(q) = T q M, i.e., if the local description (F2) applies. 
Finally, if A(q) = A2(q) is of dimension one and As(q) = T q M we say that q is a tangency point and S can 
be described near q by the normal form (F3). 
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Let us take advantage of the common expressions of the normal forms (Fl), (F2), (F3), which are all of the 
type X(x, y) = (1, 0), Y(x, y) — (0, f(x, y)), in order to investigate the local behavior of g, K, and dA. 

Lemma 20 LetX(x,y) = (1,0) andY(x,y) = (0,f(x,y)) be two smooth vector fields on R 2 . Let D = {(x,y) £ 
R 2 | f(x,y) ^ 0} and g be the Riemannian metric on D having (X,Y) as an orthonormal frame. Denote by 
K the curvature of g and by dA the Riemannian density. We have 

g = dx 2 + j2 d V 2 ' 

-2(d x f) 2 + fd 2 f 
= p ■ 

dA = — dx dy. 
4.1 An example of tangency point 

A simple example of tangency point can be observed in the 2-ARS defined by 

X(x,y) = (1,0), Y(x,y) = (0,y-x 2 ). 

The pair (X, Y) appears in the classification of planar phase portraits of pairs of vector fields given by 
Davydov in [7]. For this system one has g = dx 2 + (y — x 2 )~ 2 dy 2 , and 

K _ -2 (3x 2 + y) 
(x 2 - y) 2 

The graph of K is illustrated in Figure 3. Notice that, in contrast with the behavior of the curvature in the 
Grushin plane (see Section 3.2), in this case limsup^^Q) K (q) = +oo, while we still have liminf 9 ^( ) K(q) = 
-co. 



5 The main result 

5.1 Statement 

Let M be an orientable two-dimensional manifold and let S be an orientable 2-ARS on M. Chose a positive 
oriented atlas of orthonormal frames {(^ M , X^, y' 1 )}^/ and denote by u> the volume form on iS such that 
w(A M , Y v ) = 1 on f2 M for every /i G /. As noticed in Remark 5, cu acts as a tensor on M\ Z. Define a two-form 
dA s on M \Zby the rule dA s (V(q), W(q)) = u(V, W){q). Notice that dA s = dX^ A dY^ on W \ Z for every 
At £ I. 

Fix now an orientation S of M. Recall that the choice of 3 determines uniquely a notion of integration on 
M \ Z with respect to the form dA s . More precisely, given a cL4-intcgrable function / on C M , if for every 
q € f2, S and dA s define the same orientation at q (i.e. if S(<7) = adA s (q) with a > 0), then 

[ fdA s = [ fdA s =[f \dA s \ = f f dA. 
Jn J(n,3) Jn Jn 

Let 

M ± ={qeW\Z | fj,e I,±E(X»,Y u )(q) > 0}. 

Then / n fdA s = ± J Q fdA if f2 c A/ ± . 

For every e > let M £ = {q E M \ d(q,Z) > e}, where d(-,-) is the almost-Ricmannian distance (see 
equation (4)). We say that K is S-integrable if 

lim / K dA s 

^° Jm e 

exists and is finite. In this case we denote such limit by J KdA s . 
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Figure 3: 



Theorem 21 Let M be a compact oriented two-dimensional manifold. For a generic oriented 2-ARS on M 
such that no tangency point exists, K is S-integrable and 



where \ denotes the Euler characteristic. 

Theorem 21 is proved in Section 5.2. For a generic trivializable 2-ARS without tangency points one can show, 
thanks to topological considerations (see Section 5.4), that \(M + ) = x(M~). As a consequence, we derive the 
following result. 

Corollary 22 Let M be a compact oriented two-dimensional manifold. For a generic trivializable 2-ARS on 
M without tangency points we have 



Remark 23 In the results stated above, the hypothesis that there are not tangency points seems to be essential. 
Technically, the difficulty comes when one tries to integrate the Hamiltonian system given by the Pontryagin 
Maximum Principle applied to a system written in the normal form (F3) . However it is our hope to extend the 
Gauss-Bonnet formula even in presence of tangency points, using a more general approach. 

It is anyway interesting to notice that the hypotheses of Corollary 22 are never empty, independently of M. 
Indeed: 

Lemma 24 Every compact orientable two-dimensional manifold admits a trivializable 2-ARS satisfying the 
generic conditions of Proposition 14 and having no tangency points. 

The proof of Lemma 24 is given in Section 5.5. 
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5.2 Proof of Theorem 21 



As a consequence of the compactness of M and of Lemma 17 one easily gets: 

Lemma 25 Let M be compact and oriented. For a generic 2-ARS S on M the set Z is the union of finitely 
many curves diffeomorphic to S . Moreover, there exists Eo > such that, for every < £ < Eo, the set M\M e 
is homeomorphic to Z x [0, 1]. Under the additional assumption that M contains no tangency point, Eq can be 
taken in such a way that dM s is smooth for every < e < Eo . 

Fix M and a 2-ARS S as in the statement of Theorem 21. Thus, S can be described, around each point of 
Z, by a normal form of type (F2). 

Take eo as in the statement of Lemma 25. For every e € (0, eo), let M± = M ± n M s . By definition of dA s 
and M ± , 

/ KdA s = { KdA - { KdA. 

JM C JM+ JMs 

The Gauss-Bonnet formula asserts that for every compact oriented Riemannian manifold (N, g) with smooth 
boundary dN, we have 



/ KdA + I k g ds = 2ttx(N), 

JN JdN 



where K is the curvature of (N,g), dA is the Riemannian density, k g is the geodesic curvature of dN (whose 
orientation is induced by the one of N), and ds is the length element. 

Applying the Gauss-Bonnet formula to the Riemannian manifolds (M^,g) and (M~,g) (whose boundary 
smoothness is guaranteed by Lemma 25), we have 

KdA, = 2tt(x(M+) - X (M-)) - f k g ds + f k g ds. 

Thanks again to Lemma 25, x(M ( j t ) = x(M ± ). We are left to prove that 

lim [ / kgds - / kgds) =0. (10) 



\JdM+ JdM e 

Fix q € Z and a (F2)-type local system of coordinates (x, y) in a neighborhood U q of q. We can assume 
that Uq is given, in the coordinates (x, y), by a rectangle [—a, a] x [—6, b], a, b > 0. Assume that e < a. Notice 
that ZnU q = {0}x [-6, b] and dM e nU q = {-£, e} x [-6, b]. 
We are going to prove that 

kg ds = O(e). (11) 

Then (10) follows from the compactness of Z. (Indeed, {— e} x [—b,b] and {e} x [—6,6], the horizontal edges 
of dU q , are geodesies minimizing the length from Z. Therefore, Z can be covered by a finite number of 
neighborhoods of type U q whose pairwise intersections have empty interior.) 

Without loss of generality, we can assume that M + n U q = (0, a] x [—6,6]. Therefore, M+ induces on 
8M+ = {e} x [—6, 6] a downwards orientation (see Figure 4). The curve s i— ► c(s) = (e, y(s)) satisfying 

c(s) = -F(c(s)) , c(0) = (e,0), 

is an oriented parametrization by arclength of dM^~, making a constant angle with X. Let (Oi,02) be the dual 
basis to (X,Y) on U q n A/+, i.e., 0\ = dx and 9 2 = x^e'^^dy. According to [14, Corollary 3, p. 389, 
Vol. Ill], the geodesic curvature of <9M+ at c(s) is equal to A(c(s)), where A G h}(U q ) is the unique one-form 
satisfying 

d9i = A A 2 , cZ6» 2 = -A A 6»i . 
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Figure 4: 



A trivial computation shows that 

X = d x (x- 1 e-^ x ' y) )dy. 

Thus, 

k g (c(s)) = -S^arV^M)) (d2/(F))(c( S )) = i +a x 4>(e,y(s))- 
Denote by L\ and L2 the lengths of, respectively, {e} x [0,6] and {e} x [—6,0]. Then, 

rL 2 



fcgds = / k g (c(s))ds 

dM+nU q J-Li 



L '2 / 2 



- + d x (j)(e, (s)) ) ds 



where the last equality is obtained taking y = y(—s) as new variable of integration. 

We reason similarly on dM~ OU q , on which M~ induces the upwards orientation. An orthonormal frame 
on AI~ n U q , oriented consistently with M, is given by (X,—Y), whose dual basis is (61,— 62). The same 
computations as above lead to 

rb 

-Ay. 



[ k g ds= [ ( - - d x (j>(-e,y) J — ^ 

JdM-nu q J-b\ £ J ^ 



Define 

F(e,y) = (l+ed x <f>(e,y))e-« e »\ (12) 

Then 

I k g ds- f k g ds = \ [ (F(e,y) - F(-e,y)) dy. 

JdM+nU q JdM'nUg £ J-b 

By Taylor expansion with respect to e we get 

F(e, y) - F(-e, y) = 2d e F(0, y)e + 0(e 3 ) = 0(e 3 ) 
where the last equality follows from the relation d e F(0,y) — (see equation (12)). Therefore, 



/ kgds — / k g ds = 0(e), 

JdM+nUr, JdMrnU„ 



and (11) is proved. 
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5.3 A counterexample in the non-generic case 

In this section we justify the assumption that the 2-ARS is generic, by presenting an example of 2-ARS such 
that the conclusion of Theorem 21 does not hold, although Z is smooth and A(q) is transversal to Z at every 
point q of Z. 

Let M be the two-dimensional torus [— 7r,7r] x [— tt, tt] with the standard identifications, and consider the 
trivializable 2-ARS associated with the vector fields 

X(x, y) = (1, 0), Y(x, y) = (0, 1 - cos(x)). 

In this case Z is the circle {0} x [— tt, tt] and one among M + and M~ is empty (say M~). Notice that the 
generic condition (ii) is not verified since [X, Y](q) = at every q G Z. By Lemma 20 we have 

—j—dx dy, 

1 — cos(x) 

cos(x) — 2 
2sin(f) 2 ' 

on M \ Z. Thus, / KdA s = J M \ Z KdA = -co. 

5.4 Trivializable 2-ARSs 

The aim of this section is to characterize topologically trivializable 2-ARSs having no tangency point. 

Lemma 26 Let M be orientable. For a generic trivializable 2-ARS on M without tangency points the Euler 
characteristics of M + and M~ are equal. 

Proof. Let us consider on M a notion of angle, induced by any fixed, globally defined, Riemannian metric go. 
For every 6 in S 1 , denote by Re : TM — > TM the corresponding rotation of angle 6. 

Since, by hypothesis, the map Z 3 q t—> A(g) is a one-dimensional distribution everywhere transversal to 
the smooth submanifold Z, then we can define a smooth function <j) : Z — > (0, tt) such that R ( p^(A(q)) = T q Z 
for every q g Z. Let 6 : M X [0, 1] — > S 1 be a C°° map such that 

#|mx{o} = 
8\zx{i} = <i>- 

For every t £ [0, 1] define a vector field X t on M by the relation 

Xt{q) = Re {q ,t)X{q). 

Then X\ is a smooth vector field tangent to Z at every point of Z. Notice that {q £ M | X\{q) = 0} = {q £ 
M | 1(5) = 0} C Z. Moreover, under the generic assumption that the zeroes of X are non-degenerate, the 
same is true for those of X\ . 

Consider now the manifold M + obtained by gluing smoothly two copies of M + along Z. Since X\ is tangent 
to Z y the vector field X\, obtained as a double copy of Xi\ M +, is well defined, continuous, and has isolated 
zeroes. Thus, is equal to the sum of the indices of the zeroes of X\ . Notice that X\ (q) = if and only if 

q G Z and Xi(q) = 0. Moreover, the index of Ai at q is equal to that of X\, since the latter is non-degenerate. 
The same reasoning on M~ and X\\ M - shows that the Euler characteristic of M~, obtained by gluing two 
copies of M~ along Z, is again equal to the sum of the indices of the zeroes of X\, i.e., to x(M). Therefore, 

,(M+) = ^ = ^ = ^= X (M-). . 



dA = 

K = 
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Figure 5: 



5.5 Construction of trivializable 2-ARSs with no tangency points 

In this section we prove Lemma 24, by showing how to construct a trivializable 2-ARS with no tangency points 
on every compact orientable two-dimensional manifold. 

For the torus, an example of such structure is provided by the standard Riemannian one. The case of a 
connected sum of two tori can be treated by gluing together two copies of the pair of vector fields X and Y 
represented in Figure 5A, which arc defined on a torus with a hole cut out. In the figure the torus is represented 
as a square with the standard identifications on the boundary. The vector fields X and Y are parallel on the 
boundary of the disk which has been cut out. Each vector field has exactly two zeros and the distribution 
spanned by X and Y is transversal to the singular locus. Examples on the connected sum of three or more tori 
can be constructed similarly by induction. The resulting singular locus is represented in Figure 5B. 

We are left to check the existence of a trivializable 2-ARS with no tangency points on a sphere. A simple 
example can be found in the literature and arises from a model of control of quantum systems (see [5, 6]). 
Let M be a sphere in R 3 centered at the origin and take X(x,y,z) = (y, — x, 0), Y(x,y,z) = (0,z,—y) as 
orthonormal frame. Then X (respectively, Y) is an infinitesimal rotation around the third (respectively, first) 
axis. The singular locus is therefore given by the intersection of the sphere with the plane {y = 0} and none of 
its points is tangency (see Figure 6). Notice that the generic conditions given in Proposition 14 are satisfied. 



References 

[1] A. A. Agrachev, A "Gauss-Bonnet formula" for contact sub-Ricmannian manifolds, Dokl. Akad. Nauk, 
381 (2001), pp. 583-585. 



18 




[2] A. A. Agrachev, Yu.L. Sachkov, Control Theory from the Geometric Viewpoint, Encyclopedia of 
Mathematical Sciences, 87, Springer, 2004. 

[3] A. Bellaiche, The tangent space in sub-Ricmannian geometry, in Sub-Riemannian geometry, edited by 
A. Bclla'iche and J.-J. Risler, pp. 1-78, Progr. Math., 144, Birkhauser, Basel, 1996. 

[4] U. Boscain, B. PlCCOLi, A short introduction to optimal control, in Controle non lineaire et applications, 
edited by T. Sari, pp. 19-66, Travaux en cours, Hermann, Paris, 2005. 

[5] U. Boscain, G. Charlot, J. -P. Gauthier, S. Guerin, H.-R. Jauslin, Optimal control in laser- 
induced population transfer for two- and thrcc-lcvcl quantum systems, J. Math. Phys., 43 (2002), pp. 2107- 
2132. 

[6] U. Boscain, T. Chambrion, G. Charlot, Nonisotropic 3-level quantum systems: complete solutions 
for minimum time and minimal energy, Discrete Contin. Dyn. Syst. Ser. B, 5 (2005), pp. 957-990. 

[7] A. A. Davydov, Qualitative Theory of Control Systems, Transl. Math. Monogr., Amer. Math. Soc, Prov- 
idence, RI, 1994. 

[8] B. Franchi, E. Lanconelli, Une metrique associee a une classe d'operateurs elliptiques degeneres, Con- 
ference on linear partial and pseudodiffercntial operators (Torino, 1982), Rend. Sem. Mat. Univ. Politec. 
Torino, Special Issue (1984), pp. 105-114. 

[9] B. Franchi, E. Lanconelli, Holder regularity theorem for a class of linear nonuniformly elliptic oper- 
ators with measurable coefficients, Ann. Scuola Norm. Sup. Pisa CI. Sci., 10 (1983), pp. 523-541. 

[10] V.V. Grusin, A certain class of hypoelliptic operators (Russian), Mat. Sb. (N.S.), 83 (125) 1970, pp. 456- 
473. English translation: Math. USSR-Sb., 12 (1970), pp. 458-476. 

[11] V.V. Grusin, A certain class of elliptic pseudodifferential operators that are degenerate on a submanifold 
(Russian), Mat. Sb. (N.S.), 84 (126) 1971, pp. 163-195. English translation: Math. USSR-Sb., 13 (1971), 
pp. 155-185. 

[12] F. Jean, Uniform estimation of sub-Ricmannian balls, J. Dynam. Control Systems, 7 (2001), pp. 473-500. 



19 



[13] L.S. PONTRYAGIN, V.G. Boltyanskii, R.V. Gamkrelidze, E.F. Mishchenko, The Mathematical 
Theory of Optimal Processes, Interscicncc Publishers John Wiley and Sons, Inc, New York-London, 1962. 

[14] M. Spivak, A Comprehensive Introduction to Differential Geometry, Berkeley Publish or Perish, 1979. 

[15] M. Vendittelli, G. Oriolo, F. Jean, J. P. Laumond, Nonhomogcneous nilpotent approximations for 
nonholonomic systems with singularities, IEEE Trans. Automat. Control, 49 (2004), pp. 261-266. 



20 



